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ABSTRACT 

SN 1984A shows unusually large expansion velocities in lines from freshly syn- 
thesized material, relative to typical Type la Supernovae (SNe la). SN 1984A is 
an example of a group of SNe la which have very large blue-shifts of the P-Cygni 
features, but otherwise normal spectra. We have modeled several early spec- 
tra of SN 1984A with the multi-purpose NLTE model atmosphere and spectrum 
synthesis code, PHOENIX. We have used as input two delayed detonation models: 
CS15DD3 (Iwamoto et al. 1999) and DD21c (Hoflich et al. 1998). These models 
show line expansion velocities which are larger than that for a typical deflagra- 
tion model like W7, which we have previously shown to fit the spectra of normal 
SNe la quite well. We find these delayed detonation models to be reasonable ap- 
proximations to large absorption feature blue-shift SNe la, like SN 1984A. Higher 
densities of newly synthesized intermediate mass elements at higher velocities, 
v > 15000 km s -1 , are found in delayed detonation models than in deflagration 
models. We find that this increase in density at high velocities is responsible for 
the larger blue-shifts in the synthetic spectra. We show that the variations in 
line width in observed SNe la are likely due to density variations in the outer, 
high-velocity layers of their atmospheres. 

Subject headings: radiative transfer — supernovae: (SN 1984A) 

1. Introduction 

The current understanding of the progenitor system for SNe la features a sub- 
Chandrasekhar mass C+O white dwarf accreting mass from a companion star until 
reaching the Chandrasekhar mass, or by merging with another C+O white dwarf and then 
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exploding (for reviews see Livio 2000; Branch et al. 1995, and references therein). Most 
SNe la show a correlation between luminosity and light curve shape, but some SNe la 
show uncorrelated spectral properties (Branch & van den Bergh 1993; Hatano et al. 2000). 
Events like SN 1984A (Branch 1987; Barbon et al. 1989) and SN 1997bq (Garnavich et al. 
2000) show unusually large blue-shifts in the spectra (for more see Hatano et al. 2000). 
These SNe la have otherwise normal properties and the usual spectral features. They 
only appear different under close inspection. We have modeled SN 1984A using delayed 
detonation models, in which an initially sub-sonic flame (deflagration) transitions to a 
super-sonic shock (detonation). We have shown that the parameterized deflagration model, 
W7 (Nomoto et al. 1984), fits normal SNe la quite well (Lentz et al. 2000a). However, we 
found that delayed detonation models have abnormally large blue-shift velocities in their 
synthetic spectra. Thus, these delayed detonation models would seem to be ideal for SNe la 
with abnormally high expansion velocities. 



2. Computational Procedures 
2.1. PHOENIX 

We have used the multi-purpose spectrum synthesis and model atmosphere code 
PHOENIX 10.8 (see Hauschildt & Baron 1999, and references therein). PHOENIX has 
been designed to accurately include the various effects of special relativity important in 
rapidly expanding atmospheres, like supernovae. Ionization by non-thermal electrons from 
gamma-rays from the nuclear decay of 56 Ni that powers the light curves of SNe la is taken 
into account. We have updated the method for calculating the gamma-ray deposition from 
the single A-iteration technique of Nugent (1997) (based on the method of Sutherland 
& Wheeler 1984), to a solution of the spherically symmetric radiative transfer equation 
for gamma-rays with PHOENIX (Lentz 2000). We used an effective gamma-ray opacity, 
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K y = 0.06 Y e cm 2 g 1 (Colgate et al. 1980). In the models presented here we solve the NLTE 
rate equations for H I (10 levels/37 transitions), He I-II (19/37, 10/45), C I (228/1387), 
O I (36/66), Ne I (26/37), Na I (53/142), Mg II (73/340), Si II (93/436), S II (84/444), 
Ti II (204/2399), Fe II (617/13675), and Co II (255/2725). 

2.2. CS15DD3 

The model CS15DD3 was prepared from the hydrodynamical calculation of Iwamoto 
et al. (1999) by homologous expansion of the hydro model and by assuming a rise time of 
20 days after explosion to maximum light (e.g., Reiss et al. 2000; Aldering et al. 2000). A 
few layers near ~ 6000 km s^ 1 having non-monotonic velocities have been re-mapped, while 
we have included the corresponding density spike from the hydrodynamic interactions. 
The models, which extend to velocities > 100,000 km s -1 , have been truncated at 
40,000 km s" 1 for the pre-maximum spectra (13, 15, and 17 days after explosion) and for 
the post-maximum spectrum (28 days after explosion) to save model points to resolve parts 
of the atmosphere which are relevant to spectral formation. In tests, we have found that 
the excluded portions of the atmosphere do not affect the spectra. At each epoch, we have 
fit the luminosity to match the shape and color of the observations, while solving for the 
energy balance and converged NLTE rate equations. 

2.3. DD21c 

The model DD21c (Hoflich et al. 1998) has also been homologously expanded to the 
epoch of the observations with the same 20 day rise time. The hydrodynamical models have 
been extended from the largest velocity of ~ 25, 000 km s _1 to 30,000 km s" 1 using a steep 
power law and the same composition as the last model layer. Because of the steepness of the 
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power law the extra matter did not affect the spectra in tests. The extension was necessary 
to provide a better outer structure for the atmosphere to increase numerical stability. We 
have used the luminosities provided in Hoflich et al. (1998) for DD21c without modification 
or fitting, thus directly including the effects of time- dependence in the SN atmosphere. 



3. Results 

We have plotted the synthetic spectra from DD21c for days 13, 15, 17, and 28 after 
explosion against the observations from -7, -5, -3, and +8 days relative to maximum light 
respectively in Figure 1. We use the -7 day spectrum from Wegner & McMahan (1987). 
The epoch of maximum light and remaining spectra are from Barbon et al. (1989). The 
spectra fit remarkably well. The excellent agreement of the synthetic spectral color to 
the observations shows both that the calculated bolometric luminosities in Hoflich et al. 
(1998) are accurate, and that PHOENIX does an excellent job of reproducing time-dependent 
calculations using the externally calculated bolometric luminosities. The Si II feature near 
6000 A, the defining characteristic of SNe la, fits very well in both velocity (wavelength) 
range and shape. The general fit of the 28 day synthetic spectrum over 4500-5600 A is 
quite good, especially since the model is making the transition between optically thick and 
optically thin and therefore is more difficult to model. The S II feature at 5100-5400 A is 
shifted a little to the red in the pre-maximum spectra, indicating that the sulfur in the 
model should be moving faster to fit SN 1984A. The 4400 A and 5000 A Fe II features in 
the pre-maximum spectra are both too weak, but seem to have the correct velocity. The 
Ca II H+K feature in the day 13 spectrum is both a bit too weak and too slow. This 
indicates that higher abundance of calcium may be needed at these extreme velocities. The 
Ca II H+K feature is also quite sensitive to the model temperature. The velocity extent of 
the model is large enough to form the observed Ca II H+K line blue-shift. 
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Our best fits, adjusting the total bolometric luminosity to fit the color, to SN 1984A 
with CS15DD3 are plotted in Figure 2. Again the fit to the Si II feature at 6000 A is 
excellent. The S II in the pre-maximum spectra shows improvement over DD21c, but is a 
little slow on day 13 and weak on day 17. The 5000 A Fe II feature fits much better in all 
synthetic spectra than DD21c, especially in showing the broad, deep-bottomed shape of the 
feature. The pre-maximum 4400 A feature again is a bit weak and possibly a bit slow on 
day 13. The Ca II H+K feature shows the same weak and slow line as DD21c. The day 
28 spectrum fits fairly well except for the depth of the 6100 A Si II feature and the small 
continuum deficit just blueward of that feature. This model has a continuum optical depth, 
7"50oo = 1-7, which is quite low. 

Comparing our results to those of Figure 8 in Hoflich et al. (1998) (HWT98), the 
general features are the same, but the velocities of the lineshapes in HWT98 appear 
significantly slower than in our synthetic spectra. Figure 3 displays the 17 day spectrum in 
a format that facilitates direct comparison with Figure 8 of HWT98. There are differences 
in the Ca II IR triplet, which is significantly stronger in HWT98, also the silicon ratio 
(Nugent et al. 1995) is much larger in our calculations. We also find significantly more flux 
near 4000 A than do HWT98. 

To demonstrate the similarities between the two delayed detonation models and the 
contrasts to the deflagration model, W7, we have plotted the day 15 spectra for DD21c, 
CS15DD3, and W7 with the -5 day SN 1984A spectrum in Figure 4. We have also plotted 
W7 where we have artificially enhanced the metallicity of the unburned C+O layer by a 
factor of 10 (see Lentz et al. 2000b, where we discussed how progenitor metallicity variation 
could affect the silicon line), clearly even an extreme value for the metallicity can not be 
the cause of the high velocities of SN 1984A. We can see that although CS15DD3 is a bit 
better in the Fe II, S II, and 5700 A Si II features, the two delayed detonation models are 
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otherwise quite similar. The synthetic W7 spectrum, which has been carefully fit to the 
ordinary SN la SN 1994D (Lentz et al. 2000a), has the right color, but the spectral features 
are all too slow. To understand these differences, we have plotted the integrated column 
density for all three models in Figure 5. The delayed detonation models show larger column 
density for v > 15, 000 km s _1 , but are like W7 at slower velocities. In W7, the high velocity 
region is unburned C+O, but in the delayed detonation models it also includes some freshly 
synthesized intermediate mass elements including calcium, sulfur, and silicon. It is this 
increased density of intermediate mass elements at high velocities to which we attribute 
the high blue-shift features in SNe la like SN 1984A. We suggest that the variations in 
fio(Si II), the blue-shift velocity of the main Si II feature 10 days after maximum light, 
noted in Hatano et al. (2000) may be due to higher densities in the outermost layers of the 
supernova. 

Figures 6 and 7 compare the silicon and iron number densities of the three models. 
Clearly, DD21c and CS15DD3 fit the silicon feature well because they both have significant 
silicon (nearly the same) up to v — 30,000 km s -1 . The smoothly rising iron density of 
CS15DD3 produces a better fit to the Fe II features than does the more complicated profile 
of DD21c. 



4. Discussion 

We have calculated synthetic spectra from the delayed detonation models, DD21c 
and CS15DD3, to match the observed spectra of SN 1984A. The synthetic spectra have 
successfully reproduced the high blue-shift absorption features, and for DD21c, the colors 
using the provided bolometric luminosities. The extra column density in the intermediate 
mass elements causes the lines to form with larger blue-shifts. We have previously shown 
that line formation in this region of W7 is important in a study of C+O metallicity (Lentz 
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et al. 2000b). The variations in line velocities in SNe la with otherwise normal spectra 
are most likely due to the addition of a small amount of mass at high velocity, possibly 
containing freshly synthesized material. While the delayed detonation models we have used 
do show the higher line blue-shifts and high-velocity densities, and the deflagration model 
does not, we do not consider this as conclusive evidence that the variation is related to the 
different explosion mechanism. It may reflect a continuous variation in the ejecta densities 
at high velocity within a single family of models. A delayed detonation model with less 
density at high velocities, similar in outer density structure to W7, may fit the normal 
velocity SNe la with the same or better quality than W7. In any case it is clear that there 
must exist a family of "normal" SNe la and that spectral modeling is the way to understand 
them. 
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Fig. 1. — Synthetic spectra for DD21c (dot-dashed lines) plotted against observed spectra 
for SN 1984A (solid lines). 
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Fig. 2. — Synthetic spectra for CS15DD3 (dashed lines) plotted against observed spectra for 
SN 1984A (solid lines). 
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Fig. 3. — Synthetic spectra for CS15DD3 (dashed lines) plotted on the same scale of Figure 8 
of HWT98. 
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Fig. 4. — Day 15 synthetic spectra for CS15DD3 (short-dashed line), W7 (thick long-dashed 
line), W7 with the metallicity in the unburned C+O layer enhanced by a factor of 10 (Lentz 
et al. 2000b) (thin long-dashed line), and DD21c (dot-dashed line) plotted against observed 
-5 day spectrum for SN 1984A (solid lines). Clearly, metallicity effects cannot produce the 
extremely high velocities seen here. 



- 15 - 




Velocity (km/s) 



Fig. 5. — Column density for CS15DD3 (triple dot-dashed line), W7 (solid line), and 
DD21c (dot-dashed line) plotted against velocity. 
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Fig. 6. — Relative silicon number density for CS15DD3 (triple dot-dashed line), W7 (solid 
line), and DD21c (dot-dashed line) plotted against velocity. 
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Fig. 7. — Relative iron number density for CS15DD3 (triple dot-dashed line), W7 (solid 
line), and DD21c (dot-dashed line) plotted against velocity. 



